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We investigate the finite-temperature excitation spectrum in the gluon sector of SU(3) 



pure gauge theory through measurements of screening masses in correlations of loop op- 
erators. We develop the classification of such operators under the symmetry group of the 
'z-slice'. In the confined phase of the theory, we find that the spectrum dynamically realises 
the zero temperature symmetries. We observe a large thermal shift of the ++ glueball mass. 
In the deconfined phase, the spectrum distinguishes between operators coupling to electri- 
cally and magnetically polarised gluon fields. The former yields a screening mass equal to 
the Wilson-line screening mass; the latter, a method for the measurement of the magnetic 
mass in the high-temperature limit. 



1 Introduction 



The finite temperature deconfinement and chiral symmetry restoration transition in QCD 
and pure SU (3) gauge theory have been the subject of continuing attention. Many detailed 
studies of the screening lengths obtained from mesonic and baryonic screening correlation 
functions have been performed recently ]T|, [| . It has become clear that there is a qualitative 
difference between the hadronic excitation spectra above and below the phase transition 
temperature T c . For T < T c , the excitation spectrum resembles that at zero temperature — 
there are (nearly) massless pions, and all the well-known families of hadrons. For T > T c , 
on the other hand, quarks are liberated. As T increases beyond T c , the renormalised quark 
masses approach those computed in resummed perturbation theory M , and external currents 
with mesonic quantum numbers are almost always correlated through the exchange of a 
weakly interacting quark-antiquark pair (three weakly interacting quarks for baryons) 0. 
The exception is in the channel with pion quantum numbers. Here the inter-quark interaction 
remains quite strong even up to 2T C , but various detailed measurements show the absence 
of a pion pole H]. 

Nevertheless, there are strong spatial correlations between quarks || . A recent computa- 
tion || relates these correlations to a certain non-perturbative property of the gauge sector 
of the theory. This is the observation that spatial Wilson loops exhibit area law behaviour 
at all temperatures J7|. Non-perturbative physics has long been observed for T — > T c + in 
the pure gauge theory. This argument implies that the only non-perturbative behaviour 
observed in the fermion sector of QCD until now is related to these very same features. In 
an attempt to further study the gauge sector of the theory, we have performed extensive 
studies of certain screening correlation functions in the pure gauge theory. In the same way 
that these are used to study the excitation spectrum in the fermionic sector of QCD, we use 
them for a detailed study of the excitation spectrum in the pure gauge theory, both above 
and below T c . 

Detailed knowledge of one aspect of the pure gauge spectrum is available. In perturbation 
theory it is possible to compute the gluon propagator. From such a computation of the 
polarisation tensor, the position of the pole for the electric gluon propagator may be obtained. 
This is called the electric mass of the gluon, M e . The one- loop result has long been known 
||. For an SU(N C ) gauge group with Nf (massless) flavours of quarks, 

M 2 e = c e g 2 T\ c e = l -(N c + N } /2). (1) 
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In this paper we are interested in the case Nf = 0. It should be noted that, at non-zero 
temperatures, the appearance of an electric or magnetic mass in the gauge field propagator 
does not violate the usual Ward identities f|. Furthermore, it has been argued that the 
pole is a gauge-invariant quantity although defined through a gauge dependent correlation 
function ||. 

The electric gluon mass gives rise to the phenomenon of Debye screening of colour charges. 
This is usually seen in the screening of the colour-singlet potential between static quarks and 
the Debye screening mass can be shown to be close to M e in perturbation theory |T(|. It 
turns out to be somewhat easier to measure the screening mass in Polyakov line correlations, 
M-p. Again, in perturbation theory, it has been shown that Mp pa 2M e |TIJ. Close to 
T c , M-p has been measured through zero-momentum correlations of Polyakov loops [ |TT|| . 
Since this mass is defined through a gauge invariant correlation function, it can be given a 
non-perturbative meaning. It should be emphasised that screening masses of correlations 
between gauge invariant objects can only be connected with electric and magnetic gluon 
masses in perturbation theory. Thus screening masses are more general observables, which 
in the high-temperature and weak-coupling limit, and only in this limit, may be connected 
to the gluon electric and magnetic masses. This point has been dealt with more extensively 
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Although the one-loop computation referred to earlier shows that the magnetic part of 
the gauge field does not acquire a mass, it is possible to argue that a two-loop evaluation 
reveals quite the opposite. At finite temperature, therefore, the so-called magnetic mass, 
M m , must be of the form 

M m = c m g 2 T, (2) 

where c m is a dimensionless constant. However, all higher loop order computations contribute 
to the same order in M m , showing that this quantity cannot be perturbatively evaluated. 
Thus, the number c m above is a fully non-perturbative quantity. In this argument we have 
neglected possible terms in \ng. That these terms can be resummed into a power of g cannot 
be disproved, and indeed, has been conjectured. There were attempts long ago to measure 



M m /T in lattice simulations []12" |. We shall demonstrate that one of the screening masses 
we observe, in the high-temperature limit, goes to 2M m , and hence provides a measurement 
of M m in the same way that a measurement of Mp does for M e . In this study we have not 
been able to analyse the scaling behaviour of M m . 

Technically, our measurements are similar to zero-temperature determinations of glueball 
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masses. We construct spatial screening correlators between colour-singlet pure gauge sources 
which are loops of different shapes. These correlation functions are measured along the z- 
direction at finite temperatures. The transverse slice does not carry representations of the 
cubic group, as at T = 0, but of the group T>\ . We describe the representation theory of T>%, 
as it applies to our measurements, in section 2. The corresponding screening masses yield 
information on the physical excitation spectrum of the high temperature theory in much the 
same way as do the hadronic screening masses. 

The organisation of this paper is the following. In the next section we present the group 
theory relevant to our measurements. In section 3 we set out details of our simulations 
and the measurement procedures. The results are collected in the following section and 
conclusions presented in section 5. 

2 Symmetries of Screening Correlators 

In simulations of equilibrium field theories at finite temperatures, the rotational (Lorentz) 
symmetry of the T = theory is not realised. This is reflected in practice by the fact that 
one uses lattices of size N T x N 2 x N z , where N T < N < N z and N T = 1/aT (here a is the 
lattice spacing). Correlations between operators at different spatial separations measured 
on such an equilibrium system reveal the spectrum of screening lengths in the theory. Thus, 
one measures correlation functions of operators separated in the z-direction. As a result, it 
is immaterial whether or not N = N z ; and, in fact, it may be more convenient to choose 
N z » N. In these circumstances, it becomes necessary to identify the symmetry group 
of the continuum and lattice theories on the '^-slices', and perform the decomposition of 
operators into irreducible representations of this group. This section summarises the results 
for loop and link operators in a gauge theory. 

2.1 The Symmetry Group 

Since the correlations extend in the z direction, the relevant symmetry is of the 'z-slice'. In 
the infinite volume, continuum limit, the group action is that of the symmetries of a cylinder 
C = 0(2) x {1,<7 2 }. 0(2) is a non-Abelian group with two one-dimensional representations 
denoted Oi and O2 and two-dimensional representations labeled by the angular momentum I 
(taking integer values). The z component of the angular momentum can have the two values 
±1. The operator a z denotes reflections t — > —t about the x-y plane, and plays the role of a 
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parity. We denote its eigenvalues, {1, —1}, by the symbol P. The representations of 0(3), 
describing the zero temperature glueballs, reduce under C as 

i5(8)-of+s;i£ for p = (-!)' 

io(3)-0f + £i£ for P = -(-l)'. (3) 

On a finite lattice the symmetry of the z-slice is denoted by T>\ = D4 x {l,^}. T>± is the 
group of symmetries of a square and has 8 elements in 5 conjugacy classes — 

E: identity 

C4: 2 rotations about the t-axis by ±7r/2 

C\: rotation about the t-axis by n 

a: 2 reflections on the x — t and y — t planes 

a*: 2 reflections on the planes built from the t-axis and the 2 diagonals in the x — y plane. 

£>4 has four one- dimensional irreducible representations (irreps R) A\, A 2 , B>i and B 2 and 
one two-dimensional irrep E. The labelling of the irreps of the full symmetry group is as 
R PC where P stands for parity under a z and C for C-parity. The lowest few representations 
of C reduce under T>\ as 

Of - A* Of - A p 2 

l p ^E p 2 P ^ B{ + B p (4) 

3 P ^E P 4 P ^A P + A P . 

Thus, if the screening masses increase with increasing angular momentum, we find inequiv- 
alent T>^ representations only for / < 2. Higher angular momenta replicate these represen- 
tations. 

In order to compare to lattice data at T = 0, we also give the operators that transform 
according to the symmetry group Oh of a hypercubic lattice in x, y, t-space. Rotational 
symmetry makes this the spatial symmetry group |13|] . At any non-zero value of T, there is 
no symmetry argument for the operators and masses to split according to this classification. 



4 



The reduction Oh to T>\ yields 
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2.2 Representations of Loop Operators 
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Table 1: The 4-link operators 
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Table 2: The twisted 6-link operators 
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Table 3: The planar 6-link operators 

In this subsection we list the loop operators required to build the different representations 
of T)\. We give all the 4- and 6-link operators, as well as one group of 8-link operators which 
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Table 4: The bent 6-link operators 

is required to build the representation A^ + . This representation cannot be built out of the 
shorter loops. The complete list of loops used is 
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The notation U (p, v, • ■ •) denotes a product of link matrices U over a loop. The loop is 
specified from an arbitrary starting point by the directions p, v, etc., in order. The further 
definitions 
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Table 5: The combinations of 8-link operators which yield the representations shown; P 
takes the values ±1. 



T"i8± _ rp8 _i_ rpS fij\ 

fl,V, — fl,t fi,V,—fl,t fl,V,—fl,—t \' J 

turn out to be convenient. The C = 1 representations are always obtained by taking the 
real part of the loop and the C = — 1 representations by the imaginary part. 

The tables can be used to get two related pieces of information. The entries are the 
coefficients for linear combinations of loop operators required to obtain the indicated repre- 
sentation of T>4 . The combinations for 4- link operators are given in Table [IJ for twisted 6-link 
operators in Table Q for planar and bent 6-link operators in Tables [3] and ^ respectively; and 
finally for a set of 8-link operators in Table [|. For 4- and 6-link operators, the breakup of Oh 
representations (the T = symmetry group) for these sets of operators into representations 
of T>\ are also specified. For the 8-link operators this decomposition is not given. In our 
later analysis we shall require this information for the Bf + and representations built 
out of these operators. In these two cases, the containing representations of Oh are uniquely 
fixed to be E ++ and T 2 ++ . 

2.3 Link Operators and their Direct Products 

At high temperatures one may expect the observed correlations to be explained in terms of 
multiple gluon exchanges. Therefore it is necessary to perform a reduction of multigluon 
states to the irreps of D\. In order to derive the corresponding selection rules we start from 
the definition of a gluon field on the lattice — 

1 



{U^x) - UUx)} - tvlU^x) - UUx)} 
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Table 6: Irreducible components D of twofold products Di ® D 2 of representations of T)\. 
Note that the A 2 state vector cannot be combined into a colourless symmetric state and 
is therefore forbidden. 

for a gluon operator in any appropriately fixed gauge. Thus the electric gluon, A t , transforms 
according to the irrep A~ of V\ and the magnetic gluons, A x and A y , according to the irrep 
E+- of V\. 

The reduction of the direct product representations is straightforward. For two-gluon 
states the results are listed in Table |6|. Since each gluon is a colour octet, there is only one 
symmetric combination of two gluons to a colour singlet. Hence the A^ + state of Table ^| 
cannot be combined into a colourless state. However, the Af + , E~ + , Bi + , and -B^" + are 
allowed as two gluon states. The reduction of the three gluon states is listed in Table [7]. 
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Table 7: Irreducible components D of threefold products D\ <g> D 2 ® -D 3 of representations 
of V\. 

We note that the reduction of two and three gluon states in Tables || and [7] does not 
take into account the possible relative motion of the gluons. It thus applies directly only to 
gluons that are in an orbital s-state. We were able to observe an A 2 + state (see section 4), 
which does not appear in Tables |6] and 0. It is easy to see that it can be built of two gluons 
in an orbital p-state. 

The two possible channels in the Af + sector can thus be separated by splitting operators 
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into two groups — those containing links in the time direction, and those not. These two sets 
should give, in general, unequal screening masses. In the g(T) — > limit, these screening 
masses are then equal to 2M e and 2M m . In this limit the Bf + and screening masses 
should also yield 2M m \ although for any finite T there should be splittings between these 
three screening masses. Similarly, the masses should be asymptotically degenerate with 
these masses. In the T — > oo limit, obvious cross checks should be provided by the screening 
masses in the E h and the negative C-parity channels. 



3 Simulations and Measurements 

Simulations were carried out in pure £77(3) gauge theory at two temperatures. One was 
chosen to be in the confined phase, at about 0.75T C , and the other in the deconfined phase, 



at about 1.5T C . Since the phase transition for N T = 6 occurs at f3 c ~ 5.9 [13 , we chose 
to work with 8 x 16 3 and 4 x 16 3 lattices at fixed (3 = 5.93. Scaling was investigated by a 
separate run on a 6 x 16 3 lattice at (3 = 6.101. Both the N T = 4 and 6 lattices were at the 
same physical temperature in the deconfined phase. Further details of the runs are listed in 
Table EL 
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13000 


1300 


8 x 16 3 


5.93 


5000 


500 


6 x 16 3 


6.101 


6300 


630 



Table 8: Lattice sizes N T x iV 3 , Wilson coupling (5, number of measurements N m and size 
Nb of one jack-knife bin for analysis. 



The simulations were performed with an £77(2) subgroup over-relaxation (OR) mixed 
with a Cabbibo-Marinari pseudo-heatbath (HB) update Hl5| . We used 4 steps of OR to 
each of HB. The HB step used a Kennedy- Pendleton SU(2) heatbath ||16|| . The time for 
4 OR + 1 HB updates is 115 /is per link on a Cray Y-MP. Measurements were performed 
on configurations separated by 50 HB steps. The auto-correlations times away from T c are 
expected to be rather small, and each of these configurations is essentially uncorrelated with 
the others. 
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3.1 Fuzzed Operators 

Past experience, at T = 0, indicates that the signal to noise ratio for loop- loop correlations 
decays extremely fast with separation. It seems necessary to use a signal enhancement 
technique in order to measure correlation functions at distances greater than 2 or 3 lattice 
spacings. Accordingly, we used loops constructed out of fuzzed links |I7|]. This is a well- 



known technique. For each configuration, the link matrices generated are considered to be 
the level I = of fuzzing. Then a link UUi) at fuzzing level I is given in terms of those at 
the previous level, / — 1, by the relation 



, (9) 

SU(3) 



+ E U l - x (i)U l -\i + u)UjT\i + v + P) (U l -\i + 2/i)) f 

where the subscript denotes a projection onto SU(3). This projection is obtained through 
the polar decomposition of the general complex 3x3 matrix, M, on the right — M = ujUH . 
H is a Hermitian matrix, u a complex number and U the required special unitary matrix. 
We let the indices fi and v in eq. (Q) range over the x, y and t directions. In each direction 
we perform fuzzing up to level (log 2 iV — 1), where iV is the lattice size. Thus, for N T = 4 
we perform fuzzing up to level 1 for links in the t direction. 

We constructed the operators listed in section 2 out of the fuzzed links obtained by this 
procedure. Note also that the exclusion of the ^-direction links from the sum in eq. (Q) 
implies that the separation in the correlation functions measured is a well-defined quantity. 
Further, since the z links do not enter any of the operators we study, the loss in fuzzing is 
probably not very important. 

3.2 Mass Measurements 

Since different operators realising the same representation of T>\ can be correlated through 
the exchange of the same state, we have measured all elements of the full correlation matrix 

C a , p {z) = (L a (0)Lp(z))-(L a )(L p ), (10) 

where a and (3 label the loop operators in a given symmetry sector. The loop operators 
were summed over all sites on a z-slice, in order to project onto zero momentum. We can 



obtain the lowest lying excitation by a variational procedure [18 ] . This involves solving the 
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generalised eigenvalue problem 

C«j>{z x )Yp = X(z u z )C a>f 3(z )Yp. (11) 

Then the quantity 

g(z) = Y a C a>p (z)Y p , (12) 

where Y belongs to the lowest eigenvalue, projects onto the correlation function with the 
smallest mass. Local masses, m(z), were extracted by solving the equation 

9(z) cosh[m{z){N z /2-z)] 
g(z-l) "' cosh[m{z)(N z /2-{z-l))Y 1 ' 

where the left hand side was obtained from the measurements. Jack-knife estimators were 

used for the local masses and their errors. If the diagonalisation procedure exactly isolates 

the lowest mass into g{z), then m(z) should be independent of z within errors. In practise 

this does not happen; and we have to identify the lowest mass through a plateau in m(z). 

The effectiveness of the diagonalisation procedure is nevertheless seen in the extension of 

this plateau region. 

We restricted ourselves to diagonalising for z\ = \ and 2 at zq = in order to determine 
the leading eigenvector Y for each set of quantum numbers. The two different sets of eigen- 
vectors were found to be reasonably similar. It is useful to recall that this eigenvector gives 
the weight of each loop to the lowest state, and hence is an analogue of the modulus squared 
of the "wavef unction" . 

3.3 Polyakov Loops 

We have also measured the zero-momentum correlation functions between Polyakov loops 

c ww {z) = (fit(*)fi(o)>-(nt)(n). (14) 

The zero-momentum projection is achieved through the use of the walls 

&(z) = J2^( x >y> z )i £l(x,y,z) = tr Y[U (x,y,z,t). (15) 

xy t 

The conjugate wall operator is defined by replacing the link operators in eq. fllBD by 
their Hermitian conjugates. The static potential V(f) between a quark and an antiquark at 
distance f— (x, y, z) is usually measured through point to point correlation functions in the 
form 
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In perturbation theory the colour averaged potential above is expected to have the form 
exp(— M-pr)/r 2 . Expanding the exponential in eq. ([IB]) and substituting the perturbative 
form for the potential, one can then derive the appropriate form for C ww . In the long distance 
limit it is given by the exponential integral Ei(M-pz). We actually use the symmetrised form 

C ww (z) oc E x {M r z) + E 1 (M P (N Z - z)), (17) 

and determine M-p by an appropriate generalisation of a local mass. 

In the confining phase, V(r) cx rK e ff/T. The appropriate symmetrised form for the zero 
momentum correlation function is then 



C„„„(z) = A cosh 



Keff fN„ 



T V 2 



where A is some constant. It is clear that a local mass measurement yields K e ff/T and hence 
a plateau in this can be used for a measurement of the string tension at finite temperatures. 
The correlation function 

c aa (z) = <|nt(z)||n(o)|>-<|nt|><|n|>, (19) 

on the other hand, has trivial triality, and hence furnishes an alternative measurement of 
the mass in the Af + channel. 



4 Results 

The local masses obtained from eqs. ([H]) through fllTf ) are listed in Table |9|. We have 
underlined the values which we consider to be the best estimates of the asymptotic masses. 
In all our measurements an Af + mass turned out to be the lowest. The corresponding 
correlation functions were the least noisy, and the local masses could be followed out to 
distance 5 or 6. In this case it seems quite likely that we are able to identify the lowest 
mass in this channel. Apart from this, local masses at distance greater than 1 could only be 
measured with reasonable accuracy in the channels with quantum numbers A^ + , B± and 
B2 + - In these cases we could follow the local masses up to about distance 3. It is likely 
that the extraction of a mass in these latter channels is contaminated by the presence of 
higher states. There is more structure in our measurements than is seen in the group theory 
presented in section 2. This is good news, since it indicates non-trivial dynamics. 
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M 


rGp 




I 


z 
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2 


3 


4 


5 


6 


a 

O 


5 93 


A ++ 


lnral A ++ 

ixJ L ell . /l 


2 


0.87(1) 


0.78(2) 


77(4) 


7(1) 


0.46(8) 


0.5(2) 






A ++ 


local K++ 


2 


1.39(2) 


1.07(5) 


1 0(2) 


0.6(3) 










A ++ 




2 


1.23(2) 


63(3) 


58(5) 


0.51(8) 


6(1) 


0.6(2) 






A~ + 


local 


2 


1.84(2) 


1 6(1) 


1 2(5) 












B++ 


local, E ++ 


2 


1.38(2) 


1.21(6) 


1 0(1) 


0.8(3) 










B++ 

2 


local T„ ++ 


2 


1.45(2) 


1.27(6) 


1 2(1) 


1 2(4) 






A 

4: 




A ++ 


local, 1 


9 


1 012(6) 


80(3) 


0.73(5) 


8(1 ) 


7(3) 

U. 1 (O) 








A ++ 


1 n p a 1 in n t 


2 


1.36(1) 


1 27(4) 


1 3(1) 


7(2) 










AW + 

^2 


local 


2 


1.54(1) 


1.12(3) 


1.0(1) 


0.8(2) 












local 


2 


2.22(2) 


2.0(2) 














B++ 


local 


2 


2.37(2) 


1.9(2) 










6 


6.101 


At + 


local, t 


2 


0.97(1) 


0.64(4) 


0.61(8) 


0.51(13) 


0.43(11) 








A++ 


local, no t 


2 


1.063(7) 


0.96(2) 


0.97(7) 


1.1(2) 










A++ 


C a a 





1.0(1) 


0.56(1) 


0.46(2) 


0.51(4) 


0.49(7) 


0.48(6) 








local 


2 


1.40(2) 


1.00(3) 


0.83(6) 


0.8(2) 










Bt + 


local 


2 


1.82(1) 


1.6(1) 


1.3(2) 












B++ 


local 


2 


1.82(2) 


1.57(7) 


1.5(5) 









Table 9: Local masses, m(z), for the different representations tabulated as a function of 
z. We take the underlined values as our best estimate for the asymptotic masses. In the 
confined phase (N T = 8) we list the Oh representation to which the operator belongs, while in 
the deconfined phase (N T = 4 and 6) we indicate for the representation whether the operator 
contains time-like links or not. 

4.1 The At + Channel 

Note that the Af + mass obtained through the trivial-triality correlations of Polyakov loops 
is always equal to the lowest mass in the Af + channel extracted using loop operators. This 
provides a cross check on the proper extraction of this mass. 

Below T c this channel descends from the E ++ and Af representations of the Oh group. 
Grouping the operators constituting this representation of T>\ into two sets, one for each of 
the T = symmetries, we find that the lowest mass in the former group is about twice that 
in the latter. Such a splitting is also indicated by the eigenvectors, Yg, defined in eq. (|Tl"|). 
The low-temperature system thus seems to realise the T = symmetries dynamically. 

This is most easily demonstrated by the eigenvectors restricted to the two dimensional 
space spanned by plaquette correlations, where the plaquette operators could have a link 
in the Euclidean-time direction (P t = P* t + Py t ) or no links in this direction (P s = P^ y ). 
An operator in the E ++ representation of Oh can be formed by the combination — 2P S + P t . 
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We denote this in the shorthand notation (—2,1). The Af + representation of Oh is given 
by (1,1). Diagonalisation of the correlation matrix in this subspace gives the lowest mass of 
0.5 ± 0.1 with the corresponding eigenvector equal to (0.91 ± 0.03,1.00 ± 0.03). The next 
mass is roughly 1.0 and the eigenvector corresponding to this is (— 2.093±0.007,1.00±0.02). 
Similar results are obtained in the three dimensional subspace of planar 6-link operators. 

Since the T = symmetries seem to be generated dynamically, it is of interest to ask 
about the temperature dependence of glueball masses below T c . Our best estimate for the 
Af {Oh representation) mass at T = 0.75T C is 0.49(8). This mass has been measured at 
T = |9| to be 0.81(3) at (3 = 5.93 on a 24 3 x 36 lattice. Measurements at (5 = 5.90 on 12 4 
lattices verify that finite size effects near this coupling are small. Our measurement yields 

m Z(T Q = 7 5 ) TC) = °' 6 ± ( ^ + ' P = 5 - 93) - (20) 

Thus, the thermal shift is much larger than that seen in quenched hadron masses [20|. 

Above T c , the situation is quite different. There is no evidence of a splitting of the A~l + 
level according to the T = symmetry. On the other hand, a definite splitting is observed 
when the operators are grouped into two sets, one with loops containing no links in the t 
direction, and the other containing at least one such link. From the decomposition of the 
multi-gluon states (see Table ||) it is seen that these two sets correspond to the exchange of 
two magnetic and two electric gluons respectively. This interpretation is further supported 
by the agreement of the second of these two masses with Mp. It should be noted, however, 
that loops containing both space-like and time-like links can couple to magnetic and electric 
gluons with the lighter ones dominating the correlation function at large distances. For our 
parameters the electric gluons seem to be lighter and hence dominating. 

The eigenvectors also give evidence of such a decomposition. The two dimensional space 
of plaquette correlations for the N T = 6 lattice again demonstrates this very well. The 
lowest mass corresponds to the eigenvector (0.2 ± 0.2,1.01 ± 0.05), i.e., to an operator with 
large overlap with P t . The next eigenvector seems to have a larger overlap with P s than 
with P t . Similar results can also be obtained in the three dimensional space of planar 6-link 
operators. Finally, note that m/T is the same for the N T = 4 and 6 lattices, showing that 
finite lattice-spacing effects are under control. 
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4.2 Polyakov Loops 



The results for the effective string tension K e ff below T c and the mass Mp above T c are 
listed in Table [R]. For the effective string tension at a temperature T = 0.75T C we find 



K e ff(T) = 0.197(6). Although there is no measurement of \JK e ff at precisely this coupling 
at T = 0, it is possible to estimate this value from an interpolation of the data given in [|19| . 
We find 



KeffiT = 0.75T e 
\ K eff (T = 0) 



0.83(4). 



(21) 



Consequent to this we find that the ratio m(Af + )/^JK e ff is itself temperature dependent, 
changing from 3.5(2) at T = to a value of 2.5(5) at T = 0.75T C . 



N t 





I 


z 


1 


2 


3 


4 


5 


6 


7 


8 


5.93 


2 


0.071(1) 


0.043(1) 


0.041(1) 


0.040(3) 


0.039(3) 


0.036(4) 


0.043(3) 


4 


5.93 







0.400(8) 


0.47(3) 


0.51(5) 


0.66(12) 


0.5(2) 




6 


6.101 







0.241(8) 


0.24(1) 


0.34(4) 


0.34(6) 


0.33(6) 


0.3(1) 



Table 10: Local values of K e ff for N T = 8 and M-p for N T = 4 and 6. 

The simulations for N T = 4 and 6, corresponding to a temperature of about 1.5 T c , yield 
consistent values for M-p/T. From the data in the table, it is seen that we obtain 

M r = 2.0(2)T at N T = 4 and M v = 2.0(4)T at N T = 6. (22) 



These should be compared to the value 1.63(8), obtained at 2T C in |21[ from measurements of 
the static inter-quark potential, and the one- loop perturbation theory result 2 M e /T = 2g(T) 
where g(T) ~ 1. 

4.3 Other Channels 

At T = the measurement of excited glueball masses is not easy. Nor does a finite tempera- 
ture render this measurement any easier. Correlation functions in channels other than A[ + 
cannot be followed very far. As a result, evidence for a genuine plateau in the local masses is 
not conclusive. Although we are able to tentatively assign screening masses to correlations 
in these channels, more detailed work will be necessary to achieve the reliability obtained in 
the Af + channel. 
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A very interesting phenomenon is the observation of non-trivial correlations in the A 2 
sector. Asymptotically in T this mass should be given by 



where 2-n/N is the lowest non-zero mode for an electric gluon. This goes to Mp only in the 
limit of the spatial volume going to infinity (N — > oo). This formula works surprisingly well 
even at the temperature we study. For our lattice sizes (N = 16), taking M e to be half the 
lowest A~f + mass, we find m A -+ ps 0.90 for N T = 6 and 1.06 for N T = 4. 

It is also interesting that below T c , the A\ + mass originating from the E ++ representation 
of Oh seems to be degenerate with the Bf + mass originating from the same representation. 
This is consistent with other evidence that the zero temperature spectrum is realised dy- 
namically in this phase. Further investigation of these other glueball channels, at a sequence 
of temperatures, would be rewarding. 

5 Conclusion 

In this paper we have presented first measurements of screening masses for pure-glue oper- 
ators at high temperatures in SU(3) gauge theories. We have identified the symmetry of 
the 'z-slice group', T)\, and presented the reduction of loop operators under this symme- 
try. In order to understand the dynamics giving rise to the measured screening masses, we 
have also performed the reduction of these representations under the symmetries of the zero- 
temperature theory. Furthermore, we have examined the symmetries of the link operators, 
'gluons', and performed the reduction of direct products of these representations. 

Correlations in the Af + channel were the least noisy and, hence, the easiest to study. We 
were able to reach a detailed understanding of the dynamics in this channel. Our observations 
are completely consistent with the usual picture of a deconfinement phase transition in QCD. 
In the low-temperature phase one finds a (strong) perturbation of the zero temperature 
spectrum, whereas the high-temperature phase can be understood in terms of multi-gluon 
exchanges. Of course, as is well-known, at temperatures of 1.5T C , where our simulations were 
performed, the deconfined theory is not weakly interacting. 

In the confined phase the correlation functions were seen to dynamically reproduce the 
symmetries of the zero temperature theory. We observed this both in the spectrum and 
the 'wavef unctions' — the lowest mass in the Af + channel (of Pj) originated from the cor- 
relations of operators transforming under the Af + channel of the symmetry group of the 




(23) 
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zero temperature theory, and the next excited state in this channel could be identified as 
the correlations of the operators transforming as the E ++ representation. As a result it is 
possible to talk of thermal shifts in glueball masses for T < T c . We found that at T = 0.75T C 
the mass of the ++ glueball dropped to about half its zero temperature value. This is a 
larger mass shift than anything seen in the meson spectrum. 

The dynamical behaviour of these correlations in the deconfined phase are completely 
different. Both the mass-spectrum and the 'wavef unctions' indicate a splitting between 
correlations carried by the electric and magnetic polarisations. Group theoretically, the Af + 
representation of T>\ can arise from either two electric or two magnetic gluon exchanges in a 
relative s-wave state. Thus, to 0(g 2 (T)), the screening mass in this channel gives either 2M e 
or 2M m . Such a splitting is seen in the spectrum, and is correlated with the use of operators 
with and without time-like links. Furthermore, the Af + screening mass in the electric sector 
agrees very well with the measured Wilson-line screening mass. Hence, we can write 



These statements come with the usual caveat. The non-perturbative measurements of the 
two different kinds of screening masses in the Af + channel are valid at all temperatures, but 
the identification of these with twice the gluon electric and magnetic masses is valid only at 
high temperatures. 

Within the context of perturbation theory, corrections to the identifications made above 
arise from multi-gluon exchanges. Since direct products of three gluons do not contain Ai + , 
the 0(g 3 (T)) corrections can be obtained by resummation. Of course, these and higher 
corrections are not negligible at the couplings where our simulations are performed. Such 
corrections result in the mixing of electric and magnetic channels, as is seen easily when 
the perturbative diagrams are written out. Similar arguments may explain the fact, a little 
surprising at first sight, that M m turns out to be larger than M e . Their ratio is expected to 
be proportional to g{T), but we work at rather large values of g(T); it is entirely possible 
that at higher temperatures, and hence weaker couplings, one obtains a different ordering of 
these masses. 

The A^ + channel is correlated through similar two-gluon exchanges in a relative p-wave 
state. Thus, in perturbation theory at order greater than g 2 (T) the screening mass is split 
from the Af + channel. On lattices of finite spatial volumes this screening mass is split 
from the Af + even in the limit g(T) — > 0. Surprisingly, even for our simulations at 1.5T C , 




m 




(24) 
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the observed splitting between the screening masses in these two channels seems to be due 
entirely to this finite volume effect. It would be extremely interesting to continue these 
measurements to higher temperatures. 

A cross check on a perturbative interpretation would be provided by an unambiguous 
measurement of the + and screening masses. These are also correlated by the 

exchange of two magnetic gluons, and hence should be (asymptotically) degenerate with 
the larger mass measured in the A~l + channel. Present measurements do not support this 
strongly, but measurements at larger physical distances and at higher temperatures are 
clearly needed. We found the negative C-parity correlations extremely noisy. Better tech- 
niques for signal enhancement in these channels are clearly desirable, since the leading con- 
tributions to these are due to three gluon exchanges. 

In summary, we have performed measurements of screening masses in the pure gauge 
sector of QCD at finite temperatures, and identified the main change in dynamics across the 
phase transition temperature. In the low temperature phase we have measured the thermal 
shift in the ++ glueball mass. In the deconfined phase we have been able to identify screening 
masses which, in the high temperature limit, measure electric and magnetic gluon masses, 
and seen that they are not degenerate. In a sense, these define 'non-perturbative' gluon 
masses. Further studies of these screening masses are clearly of major interest. 
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